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We present a search for r sneutrino production using the Tevatron pp collision data collected with 
the CDF II detector and corresponding to an integrated luminosity of 1 fb _1 . We focus on the 
scenarios predicted by the R-parity violating (RPV) supersymmetric models in which r sneutrinos 
decay to two charged leptons of different flavor. With the data consistent with the standard model 
expectations, we set the upper limits on a(pp — > u T ) x BR(? T — ¥ efx, /it, er) and use these results to 
constrain the RPV couplings as a function of r sneutrino mass. 
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Supersymmetric (SUSY) extensions of the standard 
model (SM) are among the leading candidates for a 
theory which can solve the hierarchy problem and pro- 
vide a framework for unifying particle interactions pp. 
Gauge-invariant and renormalizable interactions intro- 
duced in the SUSY models can violate the conservation 
of baryon (B) and lepton (L) number and lead to a pro- 
ton lifetime shorter than the current experimental lim- 
its 0. This problem is usually solved by postulating 
conservation of an additional quantum number, R-parity 
R p = (_i)3(B-L)+2^ w here s is the particle spin [3J. 
However, models with R-parity-violating (RPV) interac- 
tions conserving spin and either B or L can also avoid 
direct contradiction with the proton lifetime upper lim- 
its [4]. Such models have the advantage that they natu- 
rally introduce lepton flavor violation and can generate 
non-zero neutrino masses and angles [5] consistent with 
neutrino-oscillation data [BJ. They can also explain the 
recently reported anomalous phase of the b to s transi- 
tion [7]. From an experimental standpoint, RPV interac- 
tions allow for single production of supersymmetric par- 
ticles (sparticles) in high-energy particle collisions and 
for sparticles to decay directly into SM particles only; 
this makes the lightest sparticle unstable and critically 
affects the experimental strategy of the SUSY searches. 
Due to their clean final-state signatures, processes of sin- 
gle slepton production followed by decay to a pair of SM 
charged leptons become promising search channels for R- 
parity violating SUSY particles [5]. 

In this Letter we report a search for a heavy r sneu- 
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trino v r produced in quark-antiquark dd annihilation 
and decaying via lepton-flavor-violating interactions into 
efi/fxr/er final states. The search is performed using 
data corresponding to an integrated luminosity of 1 fb _1 
collected in pp collisions at y^i = 1.96 TeV by the CDF II 
detector at the Tevatron. The results are analysed in the 
framework of the minimal supersymmetric extensions of 
the SM PP, where the RPV part of the superpotential 
relevant to the r sneutrino production and decay can be 
written as 

W RPV = X^UQjdk + ^L^ek (1) 

L and Q in Eq.Q are the SU(2) doublet superfields of 
leptons and quarks; e, u, and d are the 5/7(2) singlet 
superfields of leptons, u-type and d-type quarks; A' and 
A are the Yukawa couplings at the production and de- 
cay vertex respectively; the indices i, j, and k denote the 
fermion generations. We assume single-coupling domi- 
nance and the third super-generation to be the lightest. 
The couplings \' sn = 0.10 and A^ = 0.05, compatible 
with the current indirect limits [2j , are chosen as a bench- 
mark point. Heavy sneutrinos have been extensively 
searched for at LEP [9]. Recently, searches for heavy 
sneutrinos decaying into the e/i final state have been per- 
formed by the CDF [10, and DO collaborations [11 j. The 
results in this paper supersede [TO]. This analysis also 
represents the first search for lepton-flavor-violating de- 
cays of heavy sneutrinos into final states involving a third 
generation lepton, the r, at the Tevatron. 

CDF II is a general-purpose particle detector, de- 
scribed in detail elsewhere [12]. This measurement 
uses information from the central tracker |13j . calorime- 
ters [JJ] [T5] , and muon detectors [TO] for charged lepton 
reconstruction and identification. Reconstruction of pho- 
tons and 7T° mesons makes extensive use of the CES, the 
central shower maximum detector which is embedded at 
a depth of six radiation lengths within the electromag- 
netic calorimeter |14j . The luminosity is measured by 
a hodoscopic system of Cherenkov counters [T7]. The 
event geometry and kinematics are described using the 
azimuthal angle (f> around the beamline and the pseu- 
dorapidity r\ = — In tan | , where 9 is the polar angle 
with respect to the beamline. The transverse energy and 
momentum of the reconstructed particles are defined as: 
Et = -Esin0, pt = psin9, where E is the energy and p 
is the momentum. 

The data used in this measurement are collected using 
inclusive high-py electron and muon triggers which select 
high-py electron and muon candidates with \r)\ < 1.0. 
After event reconstruction, electron and muon candi- 
dates with pt > 20 GeV/c are identified using the 
procedures described in [JS]. In addition we use in- 
dependent measurements of the electron energy in the 
CES to improve the overall electron selection efficiency 
and identification of electron candidates radiating sig- 
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nificant energy due to the bremsstrahlung. The r lep- 
tons are identified via their hadronic decays as narrow 
calorimeter clusters associated with one or three charged 
tracks 19J. As the neutrino from the r decay escapes 
detection, the "visible" four-momentum of a t candi- 
date, p™, is reconstructed summing the four-momenta 
of charged particle tracks and neutral particle CES show- 
ers with a pion mass hypotheses. The resolution in p v T m 
is further improved by combining measurements of the 
track momenta and energies of the CES showers with 
the energy measurements in the calorimeter. A recon- 
structed r candidate is required to have the visible trans- 
verse energy, £7™, greater than 25 GeV and its most 
energetic track must have pr > 10 GeV/ c. The in- 



variant mass of its decay products, M™ = 




is required to be consistent with the r lepton decay: 
M™ < (1.8 + 0.0455 x (E# s /GeV-20)) GeV/c 2 , where 
the second term in the formula accounts for a degradation 
of the resolution in M™ at high energy. Reconstructed 
r candidates within AR — \J Arj 2 + A</> 2 < 0.4 of an 
identified electron or of a loosely-identified muon with 
large electromagnetic energy deposit in the calorimeter 
are excluded from the analysis. 

Events selected for the analysis are required to have 
two identified central (|?/| < 1) lepton candidates of dif- 
ferent flavor and opposite electric charge. The leptons 
have to be isolated: the extra energy measured within a 
cone of radius AR < 0.4 surrounding the leptons must 
be less than 10% of the lepton energy. Events with lep- 
tons consistent with a photon conversion or a cosmic ray 
hypothesis are removed from the analysis sample [18J. 

Signal and background studies are performed using 
Monte Carlo (MC) samples generated by pythia [2U] 
with the Tune A of CTEQ5L parton distribution func- 
tions |21j. The detector response is simulated with a 
GEANT3-based package [22]. The trigger, reconstruc- 
tion and identification efficiencies are measured using Z 
events as calibration samples [T5] , 

The predicted yield of signal events is calculated us- 
ing the next-to-leading order (NLO) pp — ¥ v r production 
cross section [23] . The total v T width is defined by the dd 
and klk decay modes as T~ = (3A' 311 -t-2A 2 3/c )M~ /16tt, 
where M~ is v T mass. 

There are several sources of background events that 
pass our analysis selections. We classify these contri- 
butions based on whether the lepton candidates recon- 
structed in these events originated from a "real" lepton 
(produced from a W or Z decay) or were a result of 
a hadron being misidentified as a lepton, lepton flavor 
mis-assignment or a secondary lepton inside a jet. We 
collectively refer to the lepton candidates of the second 
category as "fakes" and classify each contributing back- 
ground process into Type I, II and III according to the 
typical number of real and fake leptons reconstructed. 
Type I contains events with two real leptons and includes 



Z/j* -> tt, diboson (WW, WZ, ZZ) and U events. 
Type I is therefore called physics background. Type II 
includes events with one reconstructed fake lepton. They 
come from either (i) the W/Z/~f*+jet(s) events where one 
of the reconstructed leptons is in fact a jet misidentified 
as a lepton or (ii) Z/j* — » ee/nn events with one of the 
leptons misidentified as a lepton of a different flavor. The 
backgrounds in Type I and II are estimated using MC, 
and their expected event yields are normalized to the 
NLO cross sections [2~4Tl26| . Events with two fake lep- 
tons (Type III) are dominated by multi-jet events with 
two jets misidentified as leptons and 7+jets events; in 
the latter converted photon is not identified as 

such and gets reconstructed as an electron and a jet is 
misidentified as a fi or a r. The contribution of the pro- 
cesses in Type III is estimated using a data sample with 
two leptons of the same-charge and assuming no charge 
correlation between the two misidentified leptons. Both 
Type II and III are called fake background. 

The systematic uncertainties in this search arise from 
a number of sources. The uncertainty on the luminosity 
measurement is 6% [27] . Uncertainties on lepton identifi- 
cation efficiency are 3% for r's, 1% for electrons, and 1% 
for muons [15] ■ The jet-tor misidentification probabil- 
ity is measured with an accuracy of 15%. Uncertainties 
in the parton distribution functions (PDF) result in the 
systematic error on the predicted signal cross section, 
which varies from 4% to 20% and increases with the v r 
mass. Variations of the signal acceptance due to PDF 
uncertainties are less than 1%. 

We search for a signal from v T decays into lepton pairs 
in the distributions for dilepton invariant mass, M\\ (note 
that in the case of a hadronically decaying r, it is the r 
visible energy that is used to calculate the mass). The low 
mass region, 50 GeV/c 2 < M\\ < 110 GeV/c 2 , is used to 
validate the event selection and the background normal- 
ization. The observed and expected event yields in this 
region are in good agreement, as summarized in Table [I] 
With the normalization fixed, the backgrounds are ex- 
trapolated into the region M\\ > 100 GeV/c 2 , where the 
search is performed as a "blind" counting experiment. 
Figure [T] compares data distributions in M u to the SM 
expectations for each of the three channels. With no sta- 
tistically significant excesses observed, we use the data to 
set upper limits on a(pp — > v T ) x BR(? T — > e/x//xr/er). 



TABLE I: The observed and predicted event yields in the 
control region. Uncertainties on the predicted yields include 
both statistical and systematic contributions. 



Control Region 


50 GeV/c 


2 < Mn < 110 GeV/c 2 


Channel 


e^t 




er 


Physics Backgrounds 


100.2±6.5 


262.4±21.0 


309.6±24.7 


Fake Backgrounds 


9.4±3.3 


222.6±31.7 


577.8±37.0 


Total SM 


109.6±7.7 


485.0±40.9 


887.4±51.0 


Observed 


105 


477 


894 
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sensitivity for a selected M~ . The search results for 
M~ — 500 GeV/ c 2 are summarized in Table |nj Figure |2j 
shows the expected and observed 95% C.L. upper limits 
on a(pp — > v r ) x BR(i7 T — > e/x//ir/er) as a function of 

M- 

BR(? T 



. We also set 95% C.L. upper limits on X' 311 x 



ejj,/ fir / er) as shown in Table III 



TABLE II: Expected and observed number of events in e^i, 
fir, and er channels. The expected yields of v T events are 
calculated for M~ = 500 GeV/c 2 and RPV couplings Aau = 
0.10 and \ i3k = 0.05. 



100 200 300 400 500 600 700 


Channel 


e/i 




er 


Invariant Mass of e,u (GeV/c 2 ) 


M™ in ( GeV/c 2 ) 


440 


300 


310 



100 200 300 400 500 600 700 
Visible Mass of ex (GeV/c 2 ) 



FIG. 1: Expected SM and observed distributions in M\\ for 
efi, \it and er channels. Also shown is an expected v T signal 
for M~ = 500 GeV/c 2 and RPV couplings A311 = 0.10 and 



X i3 k = 0.05. 



0.1±0.02 0.2±0.03 
0.3±0.1 0.6±0.1 



Physics Backgrounds 0.03L0.01 
Fake Backgrounds 0.01L0.01 
Total SM background 0.04 ± 0.01 0.4± 0.1 0.9L0.1 
Expected signal 5.9 ± 0.1 2.0L0.1 2.7L0.1 
Observed 2 
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FIG. 2: The expected and observed 95% C.L. upper limits 
on a(pp — > v T ) x BR(5 T — > efi/ \it /er) as a function of M~ . 



In each of the channels (e/x, fj,r, er), the expected and 
observed upper limits are calculated using a Bayesian 
technique [35] at 95% credibility level (C.L.) as a function 
of M~ . For a given M~ , the limits are calculated by 
integrating the differential cross section da jdM\\ over the 
region M\\ > A/™ 111 , where M™ 1 " optimizes the search 



In conclusion, we have searched for production of a 
massive sneutrino decaying to e^t, fir, or er final states 
via R-parity violating interactions. We find the data 
consistent with the SM predictions and calculate the 
95% C.L. upper limits on the a(pp — > v T ) x BR(i/ r — > 
e/j,/ /iT/er) in the mass range up to 800 GeV/c 2 . Using 
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TABLE III: 95% C.L. upper limits on A^ 2 x BR(P T 
efi/fir/er). 



M~ (GeV/c*) 


e/i 




fir 




er 




300 


6 x 10" 


- 


4 x 10" 


4 


5 x 10" 


4 


400 


2 x 10" 


-4 


1 x 10 


-3 


9 x 10" 


-4 


500 


7 x 10" 


-4 


2 x 10 


-3 


3 x 10" 


•3 


600 


2 x 10" 


-3 


7 x 10" 


-3 


5 x 10" 


-3 


700 


8 x 10" 


-3 


2 x 10 


-2 


2 x 10" 


-2 



these cross section limits, we constrain X' 311 2 x BR(£v — > 
e/i/^tr/er) as a function of M~ . This analysis sets to the 
first Tevatron limits for lepton-flavor violating decays of 
heavy sneutrinos into final states involving a third gen- 
eration lepton. For the RPV couplings \' 3U — 0.10 and 
Ai3fe = 0.05 the observed 95% C.L. lower limits on v T 
mass are 558 GeV/ c 2 in the e/i channel, 441 GeV/ c 2 in 
the [it channel, and 442 GeV/ c 2 in the er channel. 
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